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One major characteristic of experts is intuitive judgment, which is an automatic process 
whereby patterns stored in memory through long-term training are recognized. Indeed, 
long-term training may influence brain structure and function. A recent study revealed 
that chess experts at rest showed differences in structure and functional connectivity 
(FC) in the head of caudate, which is associated with rapid best next-move generation. 
However, less is known about the structure and function of the brains of Baduk experts 
(BEs) compared with those of experts in other strategy games. Therefore, we performed 
voxel-based morphometry (VBM) and FC analyses in BEs to investigate structural brain 
differences and to clarify the influence of these differences on functional interactions. We 
also conducted graph theoretical analysis (GTA) to explore the topological organization 
of whole-brain functional networks. Compared to novices, BEs exhibited decreased and 
increased gray matter volume (GMV) in the amygdala and nucleus accumbens (NA), 
respectively. We also found increased FC between the amygdala and medial orbitofrontal 
cortex (mOFC) and decreased FC between the NA and medial prefrontal cortex (mPFC). 
Further GTA revealed differences in measures of the integration of the network and 
in the regional nodal characteristics of various brain regions activated during Baduk. 
This study provides evidence for structural and functional differences as well as altered 
topological organization of the whole-brain functional networks in BEs. Our findings also 
offer novel suggestions about the cognitive mechanisms behind Baduk expertise, which 
involves intuitive decision-making mediated by somatic marker circuitry and visuospatial 
processing. 

Keywords: amygdala, Baduk, head of caudate, intuitive judgment, resting-state functional connectivity, somatic 
marker hypothesis, voxel-based morphometry 



INTRODUCTION 

Board games such as chess have been studied by researchers from 
a variety of fields, such as economics (Levitt et al., 2011), com- 
puter science (Bouzy and Cazenave, 2001; Cai et al., 2010), and 
cognitive science (de Groot, 1965; Chase and Simon, 1973), 
because of the similarity between board games and real life in 
terms of the need to engage in decision-making and adaptive 
behavior to achieve specific goals under changing environmen- 
tal conditions. Cognitive science, in particular, has used board 
games to study cognitive expertise, as playing involves diverse 
cognitive functions such as attention, working memory, visuospa- 
tial processing, and decision-making (Chase and Simon, 1973; 
Gobet and Charness, 2006). Board-game players with the high- 
est level of skill, known as grand masters, are considered cognitive 
experts who develop the knowledge structures used in problem 
solving in a given domain through long periods of deliberate 
practice (Chase and Simon, 1973). Using these knowledge struc- 
tures, called chunks, templates, or schemas (Chase and Simon, 



1973; Gobet and Charness, 2006), experts can rapidly match the 
patterns they have learned and make faster and better decisions. 
Such chunk-driven unconscious automatic cognitive processes 
are often referred to as intuition, which is defined as the recog- 
nition of patterns or structures stored in long-term memory 
(Chase and Simon, 1973), and a number of researchers have 
proposed accounts of the mechanisms underlying intuitive judg- 
ment (Hodgkinson et al., 2009; Minavand chal et al., 2013), 
such as the following: dual-process theory, naturalistic decision- 
making (NDM), and somatic marker hypothesis (SMH). For 
example, the recognition-primed decision (RPD) model within 
the NDM approach focuses on the success of expert intuition 
(de Groot, 1965; Klein, 1998, 2008), as opposed to the heuristic - 
and-biases approach which adopts a skeptical attitude toward 
expert judgment (Kahneman and Klein, 2009). This shows how 
experts can make extremely rapid and favorable decisions by 
combining two processes: (i) an intuitive (automatic) process 
involving pattern matching based on past experience and (ii) 
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a deliberative (conscious) process involving mental simulation 
(or analysis) to imagine how a course of action will play out 
(Klein, 1993; Kahneman and Klein, 2009). The SMH emphasizes 
the influence of emotion-based signals (somatic states) emerging 
from the body, such as gut feelings on intuitive decision-making 
(Damasio, 1996; Dunn et al., 2006). Despite previous extensive 
studies on the mechanism behind intuitive expertise in board 
games, its neural basis remained largely enigmatic until the last 
two decades (Nichelli et al., 1994). Recent brain imaging studies 
during board-game play have resulted in renewed interest in the 
neural basis of cognitive expertise and have revealed brain regions 
associated with object recognition, such as the lateral occipital 
complex, occipitotemporal junction, (Bilalic et al., 2011a,b) and 
the fusiform cortex (FFC) (Bilalic et al., 2011c), with pattern 
recognition, such as the collateral sulcus (CoS) and retrosplenial 
cortex (RSC) (Bilalic et al., 2010, 201 lb), with recognition of rela- 
tions between objects, such as the supramarginal gyrus (SMG) 
(Bilalic et al., 2011a,b), and with intuitive best next-move gener- 
ation during chess play, such as the head of the caudate (HOC) 
(Wan et al, 2011, 2012). However, most neuroimaging studies 
with board-game experts have involved chess, even though Baduk 
differs fundamentally from chess in terms of the mental strategies 
involved. 

Baduk, as it is known in Korean (Go in Japanese and Weiqi in 
Chinese), is a popular board game in East Asia; it is played on 
a square board consisting of a pattern of 19 by 19 crossed lines. 
Whereas chess pieces have specific identities and functions, all 
Baduk pieces (called stones) have the same value and function. 
Rules of the game are very simple (http://english.Baduk.or.kr); 
two players, one playing with black stones and the other play- 
ing with white ones, alternately place a stone to capture as large 
an area as possible on the board by surrounding the opponent's 
stones. Despite its simple rules, Baduk is characterized by greater 
combinatorial complexity than chess due to the tremendous size 
of its game tree; the average branching factor (i.e., the num- 
ber of move choices available per turn) is approximately 200 in 
Baduk, whereas it is about 35 in chess (Keene and Levy, 1991). 
Additionally, unlike most other strategy games, Baduk cannot 
be won by a computer program, whereas computerized chess 
programs can beat even the world's best human player (Bouzy 
and Cazenave, 2001). Although chess and Baduk share common 
cognitive and affective processes, such as memory, attention, per- 
ception, and emotional regulation, the two games nonetheless 
differ in the following important ways. Given its larger game 
tree and heavy dependence on spatial positioning rather than 
on selecting pieces according to their roles, knowledge and pat- 
tern recognition with respect to spatial positioning may be more 
important in Baduk than in other strategy games (Gobet et al., 
2004). Recent neuroimaging studies on Baduk experts (BEs) have 
demonstrated increased activity in the occipitotemporal and pari- 
etal cortices, areas associated with visuospatial processing, such 
as integration of local features (Kourtzi et al., 2003) and spa- 
tial attention (Fink et al., 1996) respectively, while performing 
Baduk tasks (Chen et al., 2003; Ouchi et al, 2005). In addition 
to cognitive competences such as spatial processing, researchers 
have recently emphasized emotional processing in competitive 
board-game (Grabner et al., 2007) because based on evidence for 



the SMH (Bechara et al, 1994; Blakemore and Robbins, 2012), 
our performance (i.e., decision-making) is strongly affected by 
emotions. Thus, since board-game players experience a vari- 
ety of emotions while playing, an imbalance in the emotions 
can cause mistakes (DeGroot and Broekens, 2003). Accumulated 
evidence from neuroimaging and lesion studies implicates the 
amygdala (AMY), striatum, and orbitofrontal cortex in emotional 
processing (Phillips et al, 2003a,b), and suggests the ventro- 
medial prefrontal cortex (vmPFC), AMY, somatosensory cortex, 
and insula as regions of brain circuitry involved in the SMH 
(Damasio, 1996; Dunn et al, 2006). Particularly, the vmPFC is 
thought to play a role in generating somatic markers (Damasio, 
1996). Taken together, BEs may show differences in morphology 
and/or function in brain regions associated with spatial process- 
ing and emotion-based decision-making. However, until recently, 
there have not been studies investigating whether such specific 
differences exist in the brains of long-term trained BEs. 

Many neuroimaging studies about the learning- and practice- 
based superior performance of experts have provided evidence 
for cross-sectional differences and longitudinal changes in brain 
structure and function, known as neuroplasticity, in brain areas 
underlying specific skills. Such brain areas include the occipito- 
temporal cortex, which is associated with complex visual motions 
in jugglers (Draganski et al., 2004), the hippocampus, which 
is associated with spatial learning and memory in taxi drivers 
(Maguire et al., 2000; Spiers and Maguire, 2006), and the medial 
prefrontal cortex (mPFC)/medial orbitofrontal cortex (mOFC), 
which is associated with emotion regulation and self-referential 
processing in meditation experts (Jang et al, 2011; Kang et al, 
2013). In particular, recent studies have revealed that, compared 
to novices, chess experts demonstrate morphological differences 
in the HOC and its influence on functional circuits, showing 
a decrease in gray matter volume (GMV) and an increase in 
functional connectivity (FC) in this region during resting-state 
(Duan et al., 2012). However, whether such brain differences are 
specific to chess experts or extend to experts in other strategy 
games remains unclear. To address this issue, we used voxel-based 
morphometry (VBM) and resting-state functional connectivity 
(RSFC) analysis to compare BEs and novices in terms of GMV 
and to examine the effects of these morphological differences on 
functional brain connectivity at rest. 

RSFC analysis based on resting-state functional magnetic 
resonance imaging (rs-fMRI) reveals spontaneous or intrinsic 
functional connections of the brain, which are reflected in the cor- 
relation pattern of low-frequency blood-oxygen-level-dependent 
(BOLD) fluctuations between small regions of interest and all 
other brain regions (Fox et al, 2005). Recently, this approach 
has been extensively used in conjunction with graph theory to 
investigate the topological organization of brain networks (Wang 
et al., 2010). Graph theoretical analysis (GTA) of rs-fMRI enables 
visualization of the overall connectivity pattern across all brain 
regions and provides quantitative measurement of complex pat- 
terns of organization across a network, such as small-worldness, 
which measures global network connection efficiency. Using this 
approach, recent studies have reported differences in topological 
organization of the whole-brain functional network between per- 
sonality dimensions of extraversion and neuroticism (Gao et al., 
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2013), as well as between various brain diseases that involve cog- 
nitive impairments, such as Alzheimer's disease (Supekar et al., 
2008) and schizophrenia (Lynall et al, 2010), and healthy con- 
trols. However, the topological organization of the whole-brain 
functional network in cognitive experts is yet to be elucidated. 

We hypothesized that BEs would exhibit morphological differ- 
ences in brain regions underlying expertise in Baduk, particularly 
the occipitotemporal and parietal areas associated with visuospa- 
tial processing and spatial attention respectively, as well as the 
somatic marker circuitry involved in emotion-based decision- 
making, and that these morphological differences may be associ- 
ated with alterations in the functional circuits of these regions. We 
also predicted that the topological organization of their whole- 
brain functional network would be altered in the service of 
achieving the most efficient network for playing Baduk. To test 
our hypotheses, we employed VBM and RSFC, and further ana- 
lyzed the topological properties of the intrinsic brain connectivity 
network using a graph theoretical approach. We expect that this 
study will provide evidence for structural and functional brain 
differences in BEs, as well as offer additional insight into the 
nature of the varied and complex cognitive mechanisms that 
enable superior performance by BEs. 

MATERIALS AND METHODS 
PARTICIPANTS 

Seventeen BEs who had been training for 12.47 ± 1.50 
years were recruited from the Korea Baduk Association 
(http://english.Baduk.or.kr/). BEs experts were statistically 
matched for age, sex, and education level to 16 novices who 
knew the rules for playing Baduk. All subjects were right handed 
and had no history of neurological or psychiatric problems. 
The demographic characteristics of each group are presented in 
Table 1 . All procedures performed in this study were approved 
by the Institutional Review Board of Seoul National University 
Hospital. 

DATA ACQUISITION 

All image data were acquired using a 1.5-T scanner (Siemens 
Avanto, Germany). High-resolution anatomical images of the 
whole brain were acquired with Tl-weighted 3-D magnetization- 
prepared rapid-acquisition gradient-echo (MPRAGE) sequence 
[repetition time (TR)/echo time (TE) = 1 160/4.76 ms, flip 
angle = 15°, field of view (FOV) = 230 mm, matrix size = 256 x 
256]. rs-fMRI data were obtained via a gradient echo-planar 



Table 1 | Demographic characteristics of experts and novices. 





Experts [n = 17) 


Novices (n = 16) 


p-value 


Sex (male/female)* 


14/3 


12/4 


0.606 


Age (years)* 


17.177 (1.131) 


16.938 (1.124) 


0.547 


Education (years)* 


9.647 (2.805) 


10.688 (1.302) 


0.186 


IQ* 


93.118 (10.093) 


100.750 (12.503) 


0.062 


Training duration (years) 


12.471 (1.505) 







IQ, intelligence quotient. Values are presented as mean (standard deviation). f x 2 
test was used. * Independent t-test was used. 



imaging pulse sequence (TR/TE = 2340/52 ms, flip angle = 90°, 
FOV = 220mm, voxel size = 3.44 x 3.44 x 5mm 3 ), during 
which subjects were instructed to relax with their eyes closed 
without falling asleep. rs-fMRI scans were part of fMRI ses- 
sions, during which participants performed working memory 
tasks. Resting-state runs were performed for 4.68-min (120 vol- 
umes) prior to administration of the working memory tasks. 
Other image parameters (task-related fMRI and DTI) that are 
not related to the present study are not described herein. Based 
on visual inspection, a neuroradiologist (CHC) judged all scans 
to be excellent, without obvious motion artifacts, signal loss, or 
gross pathology. 

VOXEL-BASED MORPHOMETRY ANALYSIS 

Tl data were processed using VBM8 toolbox (http://dbm. 
neuro.uni-jena.de/vbm.html) implemented in SPM8 (http:// 
www.fil.ion.ucl.ac.uk/spm), with default parameters incorpo- 
rating the DARTEL toolbox to produce a high-dimensional 
normalization protocol (Ashburner, 2007). Images were cor- 
rected for bias-field inhomogeneities, tissue-classified into gray 
matter (GM), white matter (WM), and cerebrospinal fluid 
(CSF) based on unified segmentation from SPM8, and spatially 
normalized to the MNI space using linear (12-parameter affine) 
and non-linear transformations (warping). The nonlinear trans- 
formation parameters were calculated via the DARTEL algorithm 
(Ashburner, 2007) with an existing standard template in VBM8. 
The warped GM segments were modified to compensate for 
volume changes during spatial normalization by multiplying 
the intensity value in each voxel by the Jacobian determinants 
(modulated GMVs). Finally, the resulting GM images were 
smoothed with an 8-mm full-width at half maximum (FWHM) 
isotropic Gaussian kernel. Voxel-wise comparisons of GMV in 
the two groups were performed using two-sample f-tests. Total 
intracranial volume (TIV) was modeled as a covariate of no 
interest. TIV was calculated by summing the raw volumes of GM, 
WM, and CSF, in which each tissue volume was automatically 
generated as a text file for each subject (*_seg8.txt) in VBM8 
processing. The statistical significance of group differences was 
set at p < 0.05 using AlphaSim correction (with a combination 
of a threshold of p < 0.005 and a minimum cluster size of 340 
voxels) (Cox, 1996). Based on previous research (Duan et al., 
2012), a looser p-threshold was chosen (p < 0.005 and expected 
voxels per cluster k > 133) to detect the presence of group 
differences in the HOC. To investigate associations between 
the GMV of BEs and training duration, we employed SPM8 to 
perform voxel-wise correlation analysis between these two values 
using a multiple regression model. 

FUNCTIONAL CONNECTIVITY ANALYSIS 

rs-fMRI data preprocessing was performed using SPM8 and REST 
VI. 7 toolkit (http://www.restfmri.net/; Song et al., 2011). The 
preprocessing procedures for rs-fMRI data were performed as 
follows. After discarding the first four volumes to allow for sta- 
bilization of the BOLD signal, each subject's rs-fMRI data were 
(i) corrected for slice-timing differences, (ii) realigned to their 
first scan to correct for movement, (iii) spatially normalized to 
the MNI echo-planar imaging template in SPM8 (voxels were 
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resampled to 3 x 3 x 3 mm 3 ), (iv) spatially smoothed with a 
6-mm FWHM Gaussian kernel, (v) removed of the linear trend 
of time courses, (vi) temporally band-pass filtered (0.01-0.08 Hz), 
and (vii) conducted regression of nuisance signals (head-motion 
profiles, global signal, WM, and CSF) to correct for physiologi- 
cal noises. Regions showing significant group differences in GMV 
according to the VBM results were defined as seed regions for 
subsequent FC analysis [i.e., right AMY, right and left nucleus 
accumbens (NA); Figure 1 A] . FC maps were produced by extract- 
ing the time series averaged across voxels within each seed region 
and then computing the Pearson's correlation between that time 
series and those from all other brain voxels. Finally, correlation 
coefficients for each voxel were converted into a normal distribu- 
tion by Fischer's z transform (Fox et al., 2005). For each group, 
individual z-value maps were analyzed with a random-effect 
one-sample f-test to identify voxels with a significant positive 
correlation to the seed time series, which correlations threshold 
at p < 0.001, uncorrected, and a topological false-discovery-rate 
(FDR) correction threshold atp < 0.05 for multiple comparisons 
(Figure 2A; Chumbley et al., 2010). For between-group compar- 
isons, two-sample f-tests were used to compare z-value maps 
between experts and novices using AlphaSim correction with sig- 
nificance set at p < 0.05 (with a combination of a threshold of 
p < 0.005 and a minimum cluster size of 13 voxels for each mask 
map) (Table 2). This analysis was restricted to the voxels showing 
significant positive correlation maps for either experts or novices 
by using an explicit mask from the combined sets of the results of 
the one-sample f-tests (p < 0.05, topological FDR corrected) of 
the two groups. 

NETWORK CONSTRUCTION AND ANALYSIS 

In this study, brain networks were composed of nodes rep- 
resenting brain regions and edges representing interregional 
RSFC. To define network nodes, the Harvard-Oxford atlas 



(HOA) was employed to divide the whole brain, excluding 
the brainstem, into 110 (55 for each hemisphere) cortical and 
subcortical regions-of-interest (ROIs) (Table 3). To define the 
network edges, we calculated the Pearson correlations between 
pairs of ROIs. Correlation matrices were thresholded into binary 
networks, applying network sparsity (S) (the ratio of the number 
of existing edges divided by the maximum number of possible 
edges in a network). The sparsity threshold was normalized so 
that each group network had the same number of nodes and 
edges, allowing investigation of the relative network efficiency of 
each group (Achard and Bullmore, 2007). Given the absence of a 
gold standard for selecting a single threshold, based on previous 
studies (Wang et al, 2010; Tian et al., 2011), a continuous range 
of 0.10 < S < 0.42 with an interval of 0.01 was employed to 
threshold the correlation matrices into a set of binary matrices. 
This range of sparsity allows prominent small-world properties 
of brain networks to be observed (Watts and Strogatz, 1998); that 
is, the small-worldness of the thresholded networks was larger 
than 1.1 for all participants (Zhang et al, 2011; Gao et al., 2013). 
We calculated both global and regional network measures of 
brain networks at each sparsity threshold (Figures 3A,B, 4A,B). 
The global measures included (i) small-world parameters 
(clustering coefficient Cp, characteristic path length Lp, and 
small-worldness a) and (ii) network efficiency (local effi- 
ciency £i oc and global efficiency -E g i 0 b). The regional measures 
included three nodal centrality metrics: degree, efficiency, and 
betweenness (Rubinov and Sporns, 2010; Tian et al, 2011). 
In this study, we calculated all these metrics using GRETNA 
vl.O (https://www.nitrc.org/projects/gretna/), which is a graph- 
theoretical network analysis toolkit, with PSOM (Pipeline 
System for Octave and Matlab, (http://code.google.eom/p/psom) 
and MatlabBGL package (http://www.cs.purdue.edu/homes/ 
dgleich/packages/matlab_bgl/)- Mathematical explanations for 
each network metric are provided in the following sub-sections. 





FIGURE 1 | Brain regions showing significant group differences in gray 
matter volume and showing a correlation with Baduk training years. 

Relative to novices, Baduk experts showed significantly increased gray 
matter volume in the bilateral caudate, particularly the nucleus accumbens, 



and significantly decreased gray matter volume in the right amygdala 
(panel A). Experts showed significantly negative correlations between gray 
matter volume in the medial orbitofrontal cortex (mOFC) adjacent to the 
gyrus rectus and their training years (panels B,C). 
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A Experts 

Seed: Right amygdala 



Novices 




y = 20 y = 39 y=51 



B Experts > Novices 

Seed: Right amygdala 




Experts < Novices 

Seed: Right nucleus accumbens 




y = 51 



FIGURE 2 | Results from within- and between-group analyses of 
resting-state functional connectivity. These figures provide significantly 
positive correlation maps for the right amygdala, right nucleus accumbens, 
and left nucleus accumbens as seed regions in experts (the left column) and 



novices (the right column) (panel A). Experts displayed significantly increased 
functional connectivity between the right amygdala and left medial 
orbitofrontal cortex, and significantly decreased connectivity between the 
right nucleus accumbens and right medial prefrontal cortex (panel B). 



Table 2 | Regions showing significant group differences in gray matter volume and functional connectivity. 

Regions MNI coordinates Experts [n = 17)* Novices [n = 16) + t-value z-value 



VOXEL-BASED MORPHOMETRY RESULTS 

Experts > Novices 

Right nucleus accumbens 9,21,-6 0.428(0.027) 0.382(0.038) 3.47 3.16 

Left nucleus accumbens -11,20,-8 0.503(0.030) 0.453(0.041) 3.45 3.14 

Experts < Novices 

Right amygdala 21,2,-24 0.541(0.055) 0.573(0.031) -3.53 3.20 

Negative correlation between gray matter volumes and training duration 

Right medial orbitofrontal cortex* 3,20,-17 0.604(0.061) - 5.26 3.84 

RESTING-STATE FUNCTIONAL CONNECTIVITY RESULTS 

Experts > Novices 
Right amygdala seed 

Left medial orbitofrontal cortex -9,39,-21 0.276(0.136) 0.078(0.137) 3.68 3.32 

Experts < Novices 
Right nucleus accumbens seed 

Right medial prefrontal cortex 12,51,-6 0.102(0.120) 0.321(0.178) -3.95 3.53 

f Values are presented as mean beta values (standard deviation) for each group; mean beta values for each subject were extracted from each region using MarsBaR 
toolbox for SPIvl (http://marsbar.sourceforge.net/). 

*We performed a correlation analysis between beta value extracted from the region of and training duration of each subject using SPSS (p < 0. 001, r = —0.802,). 



Global network parameters 

To characterize the global topological organization of whole- 
brain functional network, we considered five network met- 
rics: clustering coefficient (Cp), characteristic path length (Lp), 



small-worldness (a), global efficiency (-E g i 0 b), and local effi- 
ciency (-Eioc)- Cp indicates how well neighbors of a node i 
are connected (i.e., local interconnectivity of a network). Lp 
is the shortest path length (i.e., number of edges) required to 



Frontiers in Human Neuroscience 



www.frontiersin.org 



October 2013 | Volume 7 | Article 633 | 5 



Jung et al. 



Brain differences in Baduk experts 



Table 3 | Anatomic regions-of-interest included in the network analysis 1 . 



Regions 


Abbreviation 


Classification* 


ROI index 








Left hemisphere 


Right hemisphere 


Frontal pole 


FP 


frontal 


1 


2 


Insular cortex 


Insula 


frontal 


3 


4 


Superior frontal gyrus 


SFG 


frontal 


5 


6 


Middle frontal gyrus 


MFG 


frontal 


7 


8 


Inferior frontal gyrus, pars triangularis 


IFG_PTR 


frontal 


9 


10 


Inferior frontal gyrus, pars opercularis 


IFG_POP 


frontal 


11 


12 


Precentral gyrus 


PrCG 


frontal 


13 


14 


Temporal pole 


TP 


temporal 


15 


16 


Superior temporal gyrus, anterior division 


aSTG 


temporal 


17 


18 


Superior temporal gyrus, posterior division 


pSTG 


temporal 


19 


20 


Middle temporal gyrus, anterior division 


aMTG 


temporal 


21 


22 


Middle temporal gyrus, posterior division 


pMTG 


temporal 


23 


24 


Middle temporal gyrus, temporooccipital part 


MTG TOpart 


temporal 


25 


26 


Inferior temporal gyrus, anterior division 


alTG 


temporal 


27 


28 


Inferior temporal gyrus, posterior division 


pITG 


temporal 


29 


30 


Inferior temporal gyrus, temporooccipital part 


ITG TOpart 


temporal 


31 


32 


Postcentral gyrus 


PoCG 


parietal 


33 


34 


Superior parietal lobule 


SPL 


parietal 


35 


36 


Supramarginal gyrus, anterior division 


aSMG 


parietal 


37 


38 


Supramarginal gyrus, posterior division 


pSMG 


parietal 


39 


40 


Angular gyrus 


AG 


parietal 


41 


42 


Lateral occipital cortex, superior division 


sLO 


occipital 


43 


44 


Lateral occipital cortex, inferior division 


iLO 


occipital 


45 


46 


Intracalcarine cortex 


IntraCALC 


occipital 


47 


48 


Frontal medial cortex 


FmC 


frontal 


49 


50 


Supplementary motor cortex 


SMC 


frontal 


51 


52 


Subcallosal cortex 


SubCC 


frontal 


53 


54 


Paracingulate gyrus 


ParaCG 


frontal 


55 


56 


Cingulate gyrus, anterior division 


ACG 


frontal 


57 


58 


Cingulate gyrus, posterior division 


PCG 


parietal 


59 


60 


Precuneous cortex 


PrCN 


parietal 


61 


62 


Cuneal cortex 


CN 


occipital 


63 


64 


Frontal orbital cortex 


OFC 


frontal 


65 


66 


Parahippocampal gyrus, anterior division 


aPHG 


temporal 


67 


68 


Parahippocampal gyrus, posterior division 


pPHG 


temporal 


69 


70 


Lingual gyrus 


LG 


occipital 


71 


72 


Temporal fusiform cortex, anterior division 


aTFC 


temporal 


73 


74 


Temporal fusiform cortex, posterior division 


pTFC 


temporal 


75 


76 


Temporal occipital fusiform cortex 


TOF 


temporal 


77 


78 


Occipital fusiform gyrus 


OF 


occipital 


79 


80 


Frontal operculum cortex 


FO 


frontal 


o I 


OZ 


Central opercular cortex 


CO 


frontal 


83 


84 


Parietal operculum cortex 


PO 


parietal 


85 


86 


Planum polare 


PP 


temporal 


87 


88 


Heschl's gyrus 


HG 


temporal 


89 


90 


Planum temporale 


PT 


temporal 


91 


92 


Supracalcarine cortex 


SupraCALC 


occipital 


93 


94 


Occipital pole 


OP 


occipital 


95 


96 



(Continued) 
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Table 3 | Continued 



Regions 


Abbreviation 


Classification* 


Left hemisphere 


ROI index 

Right hemisphere 


Thalamus 


Thalamus 


subcortical 


97 


98 


Caudate 


Caudate 


subcortical 


99 


100 


Putamen 


Putamen 


subcortical 


101 


102 


Pallidum 


Pallidum 


subcortical 


103 


104 


Hippocampus 


Hipp 


subcortical 


105 


106 


Amygdala 


AMY 


subcortical 


107 


108 


Nucleus accumbens 


NA 


subcortical 


109 


110 



f To define network nodes, the Harvard-Oxford atlas (HOA) was employed to divide the whole brain into 110 (55 for each hemisphere) cortical and subcortical regions 
of interest (ROIs), except the brainstem. *To facilitate data characterization and interpretation, we sorted nodes based on lobar (i.e., frontal, temporal, parietal, 
occipital, and subcortical) classification. 



A 

4.00 



Normalized C p (y) 



0.1 0.18 0.26 0.34 0.42 



Normalized E 



Loc 




Normalized L p (k) 




0.1 0.18 0.26 0.34 0.42 



Normalized E 



glob 




0.1 0.18 0.26 0.34 0.42 

FIGURE 3 | Normalized global network measures of the whole-brain 
functional network in both groups. The global measures included 
small-word parameters (panel A) and network efficiency (panel B). Asterisks 
denote significant differences (permutation-based p-value < 0.05). Significant 
differences were found in normalized path length X (permutation-based 



Small-worldness (a) 

a = y/k 



• Experts 
■ Novices 



0.1 0.18 0.26 0.34 0.42 




0.42 Normalized Lp (a) Normalized Eglob 

p-value = 0.018) and normalized global efficiency E g | 0 t, (permutation-based 
p-value = 0.008) between experts and novices (panel C). Error bars denote 
standard deviations. C P , clustering coefficient; Lp, characteristic path length; 
y, normalized clustering coefficient; X., normalized characteristic path length; 
o, small-worldness; £j 00 , local efficiency; £ g iob. global efficient. 



transfer from one node to another averaged over all pairs of 
nodes. £ g i 0 b is a measure of the capacity for parallel informa- 
tion transfer over the network, and is inversely related to Lp. 
E[ oc is a measure of the fault tolerance of the network, indi- 
cating how well each subgraph exchanges information when 
the index node is eliminated, and is related to Cp. While 
high £i oc and Cp reflect a high local specialization (called 
segregation) of information processing, high £ g i 0 b and low 



Lp express a great ability to integrate information from the 
network. 

For a given graph G with N nodes, the clustering coefficient is 
defined by Watts and Strogatz (1998) as: 



Cp(G) 



AT /-^ 



N^Dnod(l)(Dnod(!) 
i € G 



l)/2 
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Nodal Degree 



Experts 




oio ois o: 



FIGURE 4 | Regional network measures (i.e., nodal degree, nodal 
efficiency, and nodal betweenness) for experts and novices. (Panel A) 
shows values of each nodal metric over a range of thresholds in each group. 
(Panel B) shows mean values for each nodal metric across a range of thresholds 



in each group, which were superimposed on an inflated standard brain using 
BrainNet Viewer (http://www.nitrc.org/projects/bnv/). (Panel C) shows a bar plot 
of the AUC values of each nodal metric (red, frontal areas; green, temporal 
areas; blue, parietal areas; sky, occipital areas; purple, subcortical areas). 
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where -D n od(0 (see below) is the degree of a node i, and £, 
is the number of edges in G,, the subgraph consisting of the 
neighbors of a node i. The characteristic path length is defined 

by Newman (2003) as: 

L P (G) = ; j± — , 

where L« is the shortest path length between nodes i and j. To 
examine the small-world properties, the values of Cp and Lp were 
normalized as compared with those of 100 degree-matched ran- 
dom networks (y = Q eal /Q and and X = L p eal /L™ nd , a = y/X) 
before statistical analysis (Maslov and Sneppen, 2002). Typically, 
a small-world network should meet the following conditions: 
y > 1 and X « 1 (Watts and Strogatz, 1998), or 0 = y/X > 1 
(Humphries et al., 2006). The global efficiency of G is defined by 
Latora and Marchiori (2001) as: 

£elob(G) = - , 

The local efficiency of G is defined by Latora and Marchiori 
(2001) as: 

1 x - 

£loc(G) = — 2 y £glob(G,0, 

i€G 

where £ g i 0 b(G,)is the global efficiency of G;, the subgraph com- 
posed of the neighbors of a node 

Regional nodal parameters 

To investigate the regional characteristics of whole-brain func- 
tional network, we considered three nodal metrics: the nodal 
degree (-D no d)> the nodal efficiency (£ no d)> and the nodal 
betweenness (5 n od)- All these nodal metrics detect the impor- 
tance of individual nodes in the network. D no d measures the 
connectivity of a node i with all other nodes in the whole 
brain. That is, nodes with high degree interact with many other 
nodes in the network. £ noc j measures the information propaga- 
tion ability of a node i with the all other nodes in the whole 
brain. B no ^ measures the influence of a node i over informa- 
tion flow between all other nodes in the whole network. That 
is, it is the fraction of all shortest paths in the network that 
pass through a given node. The nodal degree of a node i is 
defined as: 



D„od(») = E e lp 

where ea is the i th and j th column element of the obtained 
binarized correlation matrix. The normal efficiency of a node i 
is defined as Achard and Bullmore (2007): 

s»od(o = ~ £ £ . 



The betweenness of a node i is defined as Freeman (1977): 

SnodW = 2^ g - ' 

where Syj.- is the number of shortest paths from a node; to a node 
k, and &a(t) is the number of shortest paths from a node; to a 
node k that pass through a node i within graph G. 

Statistical analysis in network parameters 

For statistical comparisons of the two groups, we calculated 
the area under the curve (AUC) for each network metric, 
which yields a summarized scalar to integrate the topologi- 
cal characteristics of brain networks over a range of thresholds 
(Figures 3C, 4C). Between-group differences in each measure 
were inferred by nonparametric permutation tests (5000 per- 
mutations) for the AUC of each global and regional measure. 
Based on a previous study (Zhang et al., 2011), we identified 
the brain regions showing significant between-group differences 
in at least one nodal metric (p < 0.05, permutation corrected). 
We also performed the Pearson correlation analyses between the 
AUC of each network metric and the duration of Baduk train- 
ing in BEs using SPSS (p < 0.05/110 for multiple comparisons 
correction). 

RESULTS 

GRAY MATTER VOLUME 

Relative to novices, BEs exhibited decreased GMV in the right 
AMY and increased GMV in the bilateral HOC, particularly the 
NA (Table 2; Figure 1 A). Significant negative correlations were 
observed between the degree of GMV reduction in the mOFC 
adjacent to the gyrus rectus and training duration in BEs (p < 
0.001, r = -0.802) (Table 2; Figures 1B,C). 

FUNCTIONAL CONNECTIVITY 

BEs showed increased FC in the right AMY seed and left mOFC, 
and decreased FC in the right NA seed and right mPFC com- 
pared to novices (Table 2; Figure 2B). We found no significant 
correlations between FC measures and training durations in BEs. 

TOPOLOGICAL ORGANIZATION OF THE WHOLE-BRAIN NETWORK 
Global network measures 

Both BEs and novices showed small-world architecture in whole- 
brain functional networks (i.e., a > 1). Compared to novices, 
BEs showed a lower normalized characteristic path length X 
and increased normalized global efficiency -E g i 0 b in the whole- 
brain functional network (Table 4; Figure 3C). No significant 
differences were found in any other global network measures. 

Regional network measures 

The groups differed with respect to nodal centrality measures 
(i.e., nodal degree, nodal efficiency, and nodal betweenness) in 
several brain regions (Table 4; Figure 5). Compared to novices, 
nodal degree in BEs showed significant increases in the right post- 
central gyrus (PocG), right inferior lateral occipital cortex (iLO), 
right thalamus, and bilateral NA, but significant decreases in the 
right superior frontal gyrus (SFG), bilateral inferior frontal gyrus 
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Table 4 | Regions showing significant differences in nodal centrality metrics between experts and novices. 

Brain regions Experts (n = 17) + Novices (n = 16) + 



p- value* 



GLOBAL NETWORK METRICS 

Experts > Novices 

Normalized global efficiency (Egi 0 t>) 
Experts < Novices 

Normalized characteristic path length (X) 
NODAL DEGREE* 
Experts > Novices 

Right postcentral gyrus 

Right lateral occipital cortex, inferior division 

Right thalamus 

Left nucleus accumbens 

Right nucleus accumbens 
Experts < Novices 

Right superior frontal gyrus 

Left inferior frontal gyrus, pars triangularis 

Right inferior frontal gyrus, pars triangularis 

Right middle temporal gyrus, posterior division 

Left lateral occipital cortex, superior division 

Right amygdala 
NODAL EFFICIENCY* 
Experts > Novices 

Right postcentral gyrus 

Right lateral occipital cortex, inferior division 

Right cingulate gyrus, posterior division 

Left thalamus 

Right thalamus 

Left nucleus accumbens 

Right nucleus accumbens 
Experts < Novices 

Right superior frontal gyrus 

Right middle temporal gyrus, posterior division 

Left lateral occipital cortex, superior division 
NODAL BETWEENNESS* 
Experts > Novices 

Right postcentral gyrus 

Right lateral occipital cortex, inferior division 

Left intracalcarine cortex 

Left parietal operculum cortex 

Right thalamus 
Experts < Novices 

Right superior frontal gyrus 

Right middle temporal gyrus, posterior division 

Right middle temporal gyrus, temporooccipital part 

Right supramarginal gyrus, posterior division 

Left temporal fusiform cortex, anterior division 

Left occipital fusiform gyrus 

Right pallidum 

Left amygdala 



0.312 (0.003) 
0.329 (0.003) 



10.216 (2.768) 
9.637 (2.401) 
9.288 (2.842) 
8.458 (2.609) 
8.012 (2.469) 

7.892 (2.966) 
7.983 (2.952) 
8.839 (2.566) 
7.092 (1.931) 
8.413 (2.554) 
7.002 (2.240) 



0.201 (0.016) 
0.198 (0.014) 
0.194 (0.024) 
0.191 (0.019) 
0.195 (0.018) 
0.190 (0.016) 
0.188 (0.015) 

0.187 (0.019) 
0.184 (0.013) 
0.191 (0.015) 



18.654 (9.765) 
17.396 (7.918) 
12.856 (6.872) 
22.417 (12.901) 
19.021 (14.113) 

13.295 (7.365) 
12.337 (3.969) 
10.113 (4.805) 
16.613 (7.972) 
7.761 (6.532) 
10.304 (5.000) 
8.257 (5.415) 
13.251 (8.945) 



0.310 (0.003) 
0.331 (0.004) 



8.404 (2.969) 
7.288 (2.087) 
7.299 (2.812) 
6.117 (2.868) 
6.376 (2.653) 

10.223 (3.530) 
9.826 (3.044) 
10.640 (3.499) 
8.706 (2.313) 
10.481 (2.537) 
8.391 (2.232) 



0.188 (0.019) 
0.184 (0.012) 
0.178 (0.030) 
0.179 (0.018) 
0.182 (0.017) 
0.171 (0.028) 
0.174 (0.026) 

0.199 (0.022) 
0.192 (0.013) 
0.201 (0.016) 



11.436 (7.248) 
11.238 (7.615) 
8.879 (5.564) 
15.055 (8.181) 
8.525 (4.291) 

18.560 (9.814) 
15.264(5.007) 
15.355 (9.793) 
22.612 (9.574) 
12.083 (7.511) 
17.177 (13.525) 
14.800 (10.731) 
20.604 (13.487) 



0.008 



0.018 



0.040 
0.002 
0.026 
0.011 
0.040 

0.025 
0.044 
0.050 
0.021 
0.013 
0.043 



0.028 
0.002 
0.038 
0.030 
0.024 
0.008 
0.029 

0.046 
0.048 
0.043 



0.011 
0.015 
0.044 
0.025 
p < 0.001 

0.043 
0.035 
0.030 
0.030 
0.045 
0.029 
0.014 
0.038 



' Values are presented as mean AUC (standard deviation) for each group. *They are p-values based on nonparametric permutation tests. 'Regions were considered 
changed in experts if they exhibited significant between-group differences (p < 0.05, permutation-corrected) in at least one nodal metric. 
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B Nodal Efficiency 




C Nodal Betweenness 




Experts > Novices ■ Experts < Novices 



FIGURE 5 | Brain regions showing significant differences in each nodal 
metric between experts and novices. Regions with significant group 
differences, including nodal degree (panel A), nodal efficiency (panel B), and 
nodal betweenness (panel C), were rendered using BrainNet Viewer 
(http://www.nitrc.org/projects/bnv/). See Table 4 for detailed information. R, 
right; L, left; Bi, bilateral; PoCG, postcentral gyrus; SFG, superior frontal 
gyurs; IFG, inferior frontal gyrus; pMTG, posterior division of middle 
temporal gyrus; sLO, superior division of lateral occipital cortex; iLO, 
interior division of lateral occipital cortex; THAL, thalamus; AMY amygdala; 
NA, nucleus accumbens; PCG, posterior cingulate gyrus; pSMG, posterior 
division of supramarginal gyrus; OF occipital fusiform gyrus; aTFC, anterior 
division of temporal fusiform cortex; PALL, pallidum; PO, parietal 
operculum cortex; TO, temporooccipital part of middle temporal gyrus; 
intraCALC, intracalcarine cortex. 



(IFG), right posterior middle temporal gyrus (pMTG), left supe- 
rior lateral occipital cortex (sLO), and right AMY (Figure 5A). 
In comparison to novices, nodal efficiency in BEs was signif- 
icantly increased in the right PocG, right iLO, right posterior 
cingulate gyrus (PCG), bilateral thalamus, and bilateral NA, but 
significantly decreased in the right SFG, right pMTG, and left 
sLO (Figure 5B). Compared to novices, nodal betweenness in BEs 
was significantly higher for the right PocG, right iLO, left intra- 
calcarine cortex, left parietal operculum cortex (PO), and right 
thalamus, while significantly lower for the right SFG, right pMTG, 
temporooccipital part of right middle temporal gyrus (TO), right 
posterior supramarginal gyrus (pSMG), left anterior temporal 
fusiform cortex (aTFC), left occipital fusiform gyrus (OF), right 
pallidum, and left AMY (Figure 5C). 



Nodal network metrics in several regions were correlated with 
training duration of BEs (p < 0.05; Figure 6), although there 
were no correlations between brain regions showing significant 
group difference in nodal network metrics and training dura- 
tions of BEs. The duration of training in BEs was positively 
correlated with nodal degree in the left SPL, but negatively cor- 
related with nodal degree in the left CN and left pTFC Nodal 
efficiency in the bilateral SPL, right anterior supramarginal gyrus 
(aSMG), and right pSMG was positively correlated with train- 
ing duration. Training duration in BEs was positively correlated 
with nodal betweenness in the left SPL, but negatively corre- 
lated with nodal betweenness in the right CN, left pMTG, and 
right PO. Figure 6 shows plots of the correlation between nodal 
metrics and training duration in BEs. However, all of these cor- 
relations did not withstand correction for multiple comparisons 
(p < 0.05/110). 

DISCUSSION 

To our knowledge, this is the first study to combine structural 
MRI and rs-fMRI to investigate the morphological differences 
in the brain of BEs, the effect of such morphological differences 
to functional circuits at rest, and the topological organization 
of the whole-brain functional network in board-game experts, 
who are treated as excellent examples of cognitive expertise. 
Four main findings emerged from this study. First, relative to 
novices, BEs showed increased GMV in the bilateral NA and 
reduced GMV in the right AMY. Additionally, the GMV in 
the mOFC was correlated with training duration. Second, BEs 
showed higher FC between the right AMY and left mOFC, and 
decreased FC between the right NA and right mPFC. Third, 
BEs showed increased global efficiency and decreased charac- 
teristic path length, implying an increase in the global integra- 
tion of the whole-brain functional network. Fourth, BEs exhib- 
ited differences in the nodal centrality metrics of many brain 
regions related to the diverse cognitive functions utilized in Baduk 
games. 

INTUITIVE EXPERTISE IN BOARD-GAME EXPERTS AND ITS 
IMPLICATION ON THE PRESENT FINDINGS 

Board games have historically been the primary focus of research 
on expertise (de Groot, 1965; Chase and Simon, 1973; Kahneman 
and Klein, 2009). Specifically, interest in how board-game experts 
make rapid and effective decisions and in the neural correlates 
associated with these intuitive judgments has been increasing in 
recent years (Duan et al., 2012; Wan et al., 2012). Such exper- 
tise or unconscious processing and linked brain changes are the 
products of prolonged experience within a domain and can- 
not be obtained through shortcuts. As it was previously stated, 
researchers have studied intuitive expertise in board-game experts 
using the NDM approach, particularly the RPD model (de Groot, 
1965; Klein, 1998, 2008; Kahneman and Klein, 2009). To describe 
how experts can make extremely rapid and good decisions, 
the RPD model combines intuitive pattern matching processes 
based on past experience and deliberative mental simulation 
processes (Klein, 2008; Kahneman and Klein, 2009). These two 
processes correspond to System 1 and System 2, respectively, 
in dual-process accounts of cognition (Hodgkinson et al., 2009; 
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FIGURE 6 | Plots of the correlations between each nodal metric and 
training duration of experts. The value of regional nodal measures 
(red, nodal degree; green, nodal efficiency; blue, nodal betweenness) 



for each BEs is shown on the x-axes and duration of Baduk training (in 
years) on the y-axes. All these correlations did not survive multiple 
comparisons correction. 



Kahneman and Klein, 2009). Interestingly, in this study, brain 
regions showing significant group difference in GMV, the NA 
and AMY, correspond to theories of the psychological and neu- 
robiological mechanisms underlying intuitive judgment. That is, 
these regions belong to the X-system (System 1), which supports 
reflexive cognitive processing in social cognition (Lieberman, 
2007) and these regions, as expected, also belong to the somatic 
marker circuitry based on the SMH. Compared with novices, BEs 
showed reduced GMV in the AMY and increased GMV in the 
NA. This finding is very interesting because of the inverse effect 
(increase/decrease) on the GMV of these two regions. Previous 
studies have reported specific correlations between regional cor- 
tical thickness and each component of cognition ( Westlye et al., 
2011), as well as a mixed pattern of increases and decreases 
in regional cortical thickness in experts (Kang et al., 2013). 
For example, meditation experts showed a thicker cortex in the 
mOFC and temporal pole, areas associated with emotional pro- 
cessing, and a thinner cortex in the parietal areas and PCG, areas 
associated with attention and self-perception when compared 
with novices (Kang et al., 2013). In addition, recent studies have 
argued for the nonlinearity of training-induced GMV changes, 
showing an initial increase followed by a decrease in regional 
GMV (Boyke et al, 2008; Driemeyer et al, 2008), and suggested 
that these changes are affected by training length and intensity 
(Takeuchi et al., 2011). Based on such previous studies, a pos- 
sible explanation for the inverse effect of the AMY and NA on 
the GMV is that each component involved in Baduk expertise 
affects a different part of the brain in a distinctive way (e.g., 
increase/decrease on GMV), or that specific processing related to 
how the amygdala shows reduced GMV may be more strongly 
involved in BEs. 



Contrary to our expectations, BEs and novices did not 
show any differences in the GMV of the occipitotemporal and 
parietal areas associated with visuospatial processing and spatial 
attention. However, we found significant group differences in 
terms of regional nodal properties of the functional brain net- 
work in these regions, particularly in the iLO, pSMG, TO, and 
pMTG, as well as in the NA and AMY, showing morpho- 
logical differences. Previous studies have reported that com- 
pared with novices, long-term chess experts (over 10 years 
of training) have morphological and functional differences in 
the HOA, showing decreased GMV and increased FC in brain 
regions known as the default-mode network (DMN) during 
resting-state in the region (Duan et al., 2012), whereas indi- 
viduals who trained for 15 weeks showed difference only in 
the activity of anatomically identified HOA but not in mor- 
phological differences in the region between before and after 
training, exhibiting increased activity during chess play (Wan 
et al., 2012). Based on previous studies and our present find- 
ings, it is, thus, conceivable that changes in brain structure, 
as well as those in brain function, through long-term train- 
ing may be required to become experts (i.e., to reach a pro- 
fessional level), and that BEs may primarily be involved in 
intuitive decision-making, which is associated with the regions 
showing both morphological and FC differences, and secon- 
darily in visuospatial processing, which is associated with the 
regions showing group differences only in the functional brain 
network. 

STRUCTURAL AND FUNCTIONAL DIFFERENCES IN THE AMYGDALA 

The observed decrease in GMV of the AMY following train- 
ing is consistent with results from previous studies that involved 
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cognitive, motor, or mental training (Boyke et al., 2008; Takeuchi 
et al., 2011; Kang et al., 2013), suggesting that a possible 
mechanism underlying such structural changes is the use- 
dependent selective elimination of synapses (Huttenlocher and 
Dabholkar, 1997). The AMY is involved in configural/holistic 
visual processing for both faces and emotional facial expressions 
(Sato et al., 201 1); it enables people to master facial affective pro- 
cessing. More specifically, the right AMY is more relevant to the 
unconscious (Morris et al, 1998) and rapid (Wright et al., 2001) 
processing of facial expressions than the left AMY. Beyond its role 
in emotional and facial processing, the AMY, in addition to the 
VMPFC, has also been proposed as one of key structures in the 
SMH and decision-making, as measured by the Iowa gambling 
task (Bechara et al., 2005), a decision-making task that requires 
implicit learning and executive function abilities. Additionally, 
animal studies have found that the AMY may participate in the 
recognition of visual patterns based on past experiences through 
direct thalamo-amygdala projection (LeDoux et al., 1989), and 
that lesions to this area may increase impulsivity in decision- 
making (Winstanley et al., 2004). Human patients with AMY 
damage have shown decreased decision-making performance 
under conditions of ambiguity and risk (Brand et al, 2007). 
Thus, GMV reduction in the AMY of BEs may be related to 
intuitive decision-making that is based on feelings (a somatic 
marker) rather than on reasoning and/or the enhanced cogni- 
tive functioning (e.g., holistic visual processing) achieved through 
long-term Baduk training. The results of our RSFC and corre- 
lation analyses support this interpretation by showing increased 
FC between the AMY and mOFC in BEs and a correlation 
between the GMV in the mOFC and training duration. Previous 
studies have demonstrated that interactions between the AMY 
and mOFC are crucial for goal-directed behavior (Holland and 
Gallagher, 2004) and emotional regulation (Banks et al., 2007). 
Lee et al. (2010) recently reported differences in the WM tract, 
the uncinate fasciculus, connecting these two regions in BEs. As 
mentioned above, the mOFC/vmPFC is an important area that 
generates somatic markers based on the SMH, and correlation 
between the GMV in the region and training duration suggests 
that morphological difference in the region may be contributed by 
Baduk training rather than pre-existing differences before train- 
ing. Contrary to this interpretation, an alternative one of the 
present finding is that this opposing pattern (reduced GMV and 
increased FC in the amygdala) may reflect a compensatory neu- 
ral mechanism pre-existing in people who later become BEs. That 
is, the observed pattern may be an endophenotype that enables 
improved somatic-marker-based (automatic) decision-making, 
and thus favorably influences (i.e., predicts) the development of 
Baduk expertise. 

STRUCTURAL AND FUNCTIONAL DIFFERENCES IN THE STRIATUM 

In contrast to the GMV in the AMY, the GMV in the NA 
was increased in BEs compared to novices. This is espe- 
cially intriguing, as a recent study with chess experts using 
the same method reported a GMV decrease in the dor- 
sal part of HOC, which is known to be involved in the 
quick generation of the best next move during chess play- 
ing (Wan et al., 2011), while we found a GMV increase in 



the ventral part of the HOC, the NA (Duan et al, 2012). 
Whereas the dorsal striatum is associated with sensorimotor 
experiences or reward-outcomes processing, the ventral stria- 
tum is associated with emotional and motivational experi- 
ences or reward-anticipation processing. Consistent with the 
present finding, Boyke et al. (2008) found increased GMV 
in the NA after juggling training. They suggested that learn- 
ing to juggle may stimulate an increase in the size of this 
region due to its role as an interface between the limbic 
and motor systems rather than any role in motor control 
(Mogenson et al, 1980). 

We also found decreased FC between the NA and mPFC 
in BEs compared with novices. Such connections are seen 
in the anterior part of the DMN (Di Martino et al, 2008; 
Duan et al., 2012; Jung et al., 2013), which is thought to 
be involved in self-referential processing and theory of mind 
(ToM) (Murdaugh et al, 2012; Reniers et al, 2012). ToM refers 
to the ability to infer the thoughts or intentions of others. 
Given the hyperconnectivity within the DMN of individuals 
at high-risk for psychosis (Shim et al, 2010), who also suf- 
fer from impaired ToM (Chung et al, 2008), a reverse pat- 
tern (decreased FC) in BEs may be due to their capacity to 
infer the opponent's intention while playing Baduk. It is also 
conceivable that the strength of this connection may reflect a 
change in their sensitivity to feedback, given that these couplings 
are thought to alter expectations in the face of negative feed- 
back (van den Bos et al., 2012). However, these ideas are just 
speculations, and additional neuroimaging studies with Baduk 
tasks are necessary to identify the physiological mechanisms 
that underlie these brain differences and to clarify the relation- 
ship between their cognitive functions and brain structure and 
function. 

TOPOLOGICAL ALTERATION IN THE WHOLE-BRAIN FUNCTIONAL 
NETWORK 

The GTA results reveal increased global efficiency and decreased 
characteristic path length in Bes compared with novices, implying 
an increase in the global integration of the whole-brain network. 
This can be reflective of effective integrity and rapid informa- 
tion propagation between and among the remote regions of the 
brain involved in the cognitive processing required for Baduk 
play in BEs (Wang et al., 2010). Thus, this finding reflects a dif- 
ference in the functional aspect of the whole-brain circuitry in 
the service of achieving the most efficient network for playing 
Baduk. 

We also found the following differences between the two 
groups in the regional nodal characteristics of many brain 
regions: increased nodal centrality metrics in nine regions, 
namely the right PocG, right iLO, right PCG, left intracal- 
carine cortex, left PO, bilateral NA, and bilateral thalamus, and 
decreased nodal centrality metrics in 12 regions, namely the right 
SFG, right pMTG, right TO of MTG, right pSMG, right pal- 
lidum, left sLO, left aTFC, left OF, bilateral IFG, and bilateral 
AMY. Interestingly, whereas the brain regions showing increased 
nodal centralities in BEs were involved primarily in implicit pro- 
cessing, such as somatic sensation (PocG and thalamus), visual 
expertise (iLO), and affective/motivational processing (NA), the 
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brain regions showing decreased nodal centralities were related 
to higher-order cognitive functions such as executive function 
(SFG, IFG, and sLO), semantic memory processing (pMTG and 
pSMG), and visual perception (TO and OF). It is speculated 
that brain regions showing differences in nodal centralities are 
important contributors in Baduk expertise, and may facilitate 
efficient exchange of information. In this context, our findings 
are consistent with the RPD model mentioned above (Klein, 
2008) Therefore, Baduk is thought to involve diverse cogni- 
tive functions with respect to both automatic and deliberative 
processes. 

Previous neuroimaging studies that focused on Baduk 
reported enhanced activation in the occipitotemporal and 
parietal cortices, during these games (Chen et al, 2003; 
Ouchi et al., 2005), as well as learning-induced differ- 
ences in activity in fronto-parietal and visual cortices (Itoh 
et al., 2008), which corresponds to our findings of differ- 
ences in the PocG and iLO. Correlation analyses between 
each nodal metric and training duration of BEs showed 
positive correlations between these two values in the pari- 
etal areas, particularly more extensive in the right than left 
hemisphere, although it did not remain significant after cor- 
rection for multiple comparison. The right SPL in terms 
of Baduk play, may contribute to spatial working memory 
(Ungerleider et al., 1998) and/or spatial attention (Fink et al., 
1996). 

Chess requires the recognition of the identity and function 
of each piece (i.e., chess-specific object and function recogni- 
tion, which is known to be associated with the occipito-parieto- 
temporal junction, OTJ; Bilalic et al., 2011a,b), whereas this 
is not needed in Baduk. This difference as a confounding fac- 
tor makes it difficult to compare or interpret the results of the 
present study with those of previous neuroimaging studies on 
chess experts. However, both Baduk and chess involve pattern 
recognition with respect to the spatial positioning of stones or 
pieces, which is an essential component in improving both games. 
Recently, Bilalic et al. (2010, 2011b) demonstrated that the CoS 
and RSC, part of the parahippocampal place area, play an impor- 
tant role in chess-specific pattern recognition. Intriguingly, we 
found increased nodal centralities in the iLO, which corresponds 
to the OTJ associated with chess-specific object and function 
recognition in studies by Bilalic et al. (2011a,b), in BEs com- 
pared to novices but not any differences in the Cos and RSC 
between the two groups. This discrepancy in results between 
the present study with BEs and previous studies with chess 
experts may stem from differences between the basic aspect of 
pattern recognition in these board games; pattern recognition 
in Baduk is only based on the shapes the stones take, while 
that in chess is based on object and function recognition of 
chess pieces and their ability to rapidly access the information 
of potential moves or move sequences for each piece. Thus, dif- 
ferences in the iLO may be linked to the recognition of shapes 
and integration of local features (Kourtzi et al., 2003) during 
Baduk play. 

Taken together, these GTA results provide such insight into 
the topological organization of the functional brain networks of 
BEs: increased functional integration across global brain regions 



and increased nodal centralities in regions associated with spatial 
attention and somatic sensation. 

LIMITATIONS 

The present study has some limitations. First, it is difficult 
to accurately and quantitatively assess the skill level of BEs. 
Although we used training duration (in years) as a proxy for 
skill level in BEs, the validity of this proxy can be challenged. 
Skill levels in BEs may be independent of training duration, 
resulting in similar skill levels for any given pair of short- and 
long-trained participants. Previous studies have described that 
the lack of any significant correlation between brain imaging 
data and the amount/intensity of training in experts may result 
from the inaccuracy of the proxy chosen in determining the 
actual extent/intensity of the individual training (Luders et al., 
2011; Kang et al, 2013). Second, as a result of the cross-sectional 
nature of this study, our findings do not allow any unambiguous 
definitive causality. Therefore, it is unclear whether the brain dif- 
ferences we found result from acquiring Baduk expertise through 
prolonged training, or if they simply reflect pre-existing differ- 
ences in brain structure and function that predict later expertise. 
Longitudinal studies will help to clarify this issue. Finally, the 
present findings with respect to FC and GTA are based on rs- 
fMRI data of BEs rather than on brain activity during game 
performance. Thus, further functional imaging studies are nec- 
essary to investigate the topological properties and FC within 
the functional brain network while the individuals actually play 
Baduk. 

CONCLUSIONS 

The current study demonstrates differences in the structure and 
the functional circuits of the AMY and NA in BEs; compared 
with novices, experts showed decreased GMV and increased FC 
with the mOFC in the AMY as well as increased GMV and 
decreased FC with the mPFC in the NA. As interfaces between 
the cognitive and affective components of the limbic cortico- 
striatal loop, the AMY and NA are involved in implicit processing 
and goal-directed adaptive behavior under changing environ- 
mental conditions. In particular, the AMY is critical for emo- 
tional and holistic visual processing, as well as emotion-based 
decision-making. Based on our hypothesis that long-term Baduk 
training would influence the structure and functional circuits 
of regions associated with the cognitive mechanisms underlying 
Baduk expertise, our findings suggest that intuitive decision- 
making, which is mediated by somatic marker circuitry such 
as the AMY and NA, is a key cognitive component of Baduk 
play. The current study also provides new evidence for differ- 
ences in the topological organization of the whole-brain net- 
work of BEs, showing increased global integration and altered 
regional nodal centralities in the regions related to visuospatial 
processing. 
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